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Possible detection of a magnetic field in T Tauri 
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Abstract. Medium-resolution {R ~ 15000) circular spectropolarimetry of T Tauri is presented. The star was 
observed twice: on November 11, 1996 and January 22, 2002. Weak circular polarization has been found in 
photospheric absorption lines, indicating a mean surface longitudinal magnetic field Sy of 160 ±40 G and 140 ±50 
G at the epoch of the first and second observations respectively. While these values are near the detection limit 
of our apparatus, we belive that they are real. In any case one can conclude from our data that By of T Tau does 
not significantly exceed 200 G, which is much less than surface magnetic field strength of the star (> 2.3 kG) 
found by Guenther et al. (1999) and Johns-KruU et al. (2000). We discuss possible reasons of this difference. 
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1. Introduction 

Classical T Tauri stars (CTTSs) are low mass pre-main se- 
quence stars which have accretion disks (Bertout, 1989). 
Magnetic fields are belived to play a crucial role in the 
evolution of CTTSs angular momentum and in the in- 
teraction between the central star and the circumstellar 
matter. Surface magnetic fields ~ 1 kG are expected from 
theoretical arguments for these objects - see e.g. Shu et al. 
(2000) and references therein. 

Meanwhile, no field has been directly detected in 
CTTSs until recently although several attempts have 
been made (Babcock 1958; Brown & Landstreet 1981; 
Johnstone & Penston 1986,1987). Recently two groups in- 
dependently succeded in detecting magnetic fields on T 
Tau, LkCa 15 (Guenther et al., 1999; Johns-KruU et al., 
2000) and BP Tau (Johns-KruU et al., 1999a,b). 

Two general methods have been used to measure the 
magnetic field of T Tau. The first one is based on the fact 
that magnetic fields produce an enhancement of line equiv- 
alent widths that vary with magnetic sensitivity (Lande 
factor and line splitting pattern). With this technique 
one can derive the product of the surface magnetic field 
strength, B, and surface filling factor /, i.e. B ■ /-value - 
see Basri et al. (1992) for details. In the case of T Tau, for 
example, Guenther et al. (1999) found B- f ^ 2.35 ±0.15 
kG, which means the average value of B itself is > 2.3 kG. 

The other direct method for detecting magnetic fields 
is to look for net circular polarization in Zeeman-scnsitivc 



lines. Generally, Zeeman a components are elliptically po- 
larized, with the components of opposite helicity split to 
either side of the nominal wavelength Aq. If there is a net 
longitudinal component of the magnetic field on the stellar 
surface, net polarization results in a shift between Zeeman- 
sensitive lines observed in right- and left-circularly polar- 
ized (RCP and LCP) light. The shift is 



AXri ~ 9.3- 10^'^ gXlB«, mA 



(1) 
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where g is the Lande factor of the transition, which we 
adopt from the VALD database (Kupka at al. 1999). Here 
the mean longitudinal magnetic field averaged over the 
star including limb darkening, B\\ , and Aq are expressed 
in kG and A respectively. 

Only upper limits of -B|| were found in the case of T 
Tau with this technique: 1000 G (Babcock, 1958) and 
810 G (Brown & Landstreet, 1981) at the 3cr level. We 
present here results of our more sensitive spectropolari- 
metric measurements of T Tau's longitudinal magnetic 
field. The following set of T Tau's atmospheric parame- 
ters were adopted: Te// ~ 5250 K, (KOV spectral type), 
log 5 = 3.73 (White & Ghez, 2001) and solar metallicity. 

2. Observations 

We observed T Tau with the 6 m telescope of the Special 
Astrophysical Observatory on November 21, 1996 and 
January 22, 2002. The Main Stellar Spectrograph (MSS - 
Panchuk, 2001) equipped with the polarimetric analyzer, 
(Najdenov & Chountonov, 1976; Chountonov, 1997) was 
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used. The analyzer splits the incoming beam into two par- 
allel beams that are offset from one another along the 
spectrograph slit by 5 arcseconds. One beam contains the 
RCP component of the spectrum and the other the LCP 
component. As a result, two spectra with opposite circu- 
lar polarization ~ 200 A in length were projected onto an 
1160 X 1024 pixel CCD. The spectral band 6300 - 6500 
A was observed in 1996 and 6600 - 6800 A in 2002. The 
width of MSS's slit was 0.5" providing a 2-pixel spectral 
resolution of i? = A/AA ~ 15000. 

Spectra were processed as follows. Dark current, night 
sky emission and detector bias as well as cosmic ray traces, 
were removed in a standard way, using routines from the 
MIDAS software package. A spectrum of a thorium-argon 
lamp was used for wavelength calibration. Then we iden- 
tifed spectral lines adopting line information from the 
VALD database (Kupka at al., 1999). The shifts between 
positions of identical lines in spectra with opposite circu- 
lar polarization in the same image, AA^;, were derived by 
means of a Fourier cross-corellation method. The method 
appears to be relatively insensitive to uncertainties in the 
continuum level position and gives the ability to evaluate 
the displacement of line profiles as a whole (Klochkova 
et al., 1996). 

We took into account systematic instrumental effects 
(SIEs) in different ways reducing the 1996 and 2002 ob- 
servations. In the first case, we used observations of the 
G9.5 III star e Tau, which presumbly has negligible field. 
Specifically, we assumed that nonzero values of AXri{X) 
in £ Tau's spectrum result from instrumental effects only 
and used them to correct AA^rvalues in the spectrum of 
T Tau - see Romanyuk et al. (1998) for more details. For 
this reason, our November 21, 1996 observations were or- 
ganized in the following way: initially we observed e Tau, 
then T Tau and then s Tau again. Approximate signal- 
to-noise ratio (S/N) of individual spectra, number of ob- 
served spectra, n, visual magnitude V, and total exposure 
time, te (in seconds), for each target are presented in Table 
1. 

Table 1. Journal of Observations 



LCP spectra is recorded in the first exposure (A^, A^) and 
in the second (A2, A2). Then we calculated 



Date 


Target 


V 


n 


te 


S/N 


22.11.1996 


e Tau 


3.5 


5 


220 


300 




T Tau 


9.9 


3 


3600 


150 


21.01.2002 


75 Tau 


5.0 


3 


1320 


300 




£ Tau 


3.5 


3 


600 


300 




T Tau 


9.9 


3 


7200 


100 


24.01.2002 


HD 30466 


7.3 


1 


1400 


200 



In the case of the 2002 observations, SIEs were re- 
moved as follows. Between exposures the phase compen- 
sator was advanced 1 /2 wave in order to reverse the sense 
of circular polarization recorded in the two spectra. For 
each target, the wavelength of a given line in the RCP and 



AA. 



A'^ - A' 



(A^ - X[) + (A^ - A^) 



which is expected to be free from the main instrumental 
effects - sec Johns-KruU ct al. (1999b) for details. 

Information on the January 2002 observations is also 
presented in Table 1, where now n stands for the number 
of pairs of individual observations of any particular star, 
between which the phase compensator was rotated. 

3. Results 

3.1. 1996 observations 

Figure 1 presents values of AAr;(A) derived from all the e 
Tau spectra observed immediately before and after obser- 
vations of T Tau. The AA^; (A) dependence was approx- 
imated with a linear function AA,.; = a ■ N + b, where 
N is the pixel number, for both observations. Wc found 



a = (-2.8 ±2.6) ■ IQ- 



per pixel, h = 10.5 ±1.4 mA 



for the open squares in Figure 1 and a = (— 5.3±3.4)-10^'^ 
mA per pixel, b = 14.4 ±1.9 mA for black triangles in 
Figure 1. 

The difference between the parameters of the fits for 
both datasets (squares and tringles) is whithin 2.5 a, i.e. 
is not significant from a statistical viewpoint. So one can 
use the resulting linear function, derived from all points, to 
remove SIEs from AA^r values in the spectra of T Tau. We 
found a = (-4.1 ±2.2) -10-3 mA per pixel, b= 12.4 ±1.2 
mA in this case - solid line in Fig.l 
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Fig. 1. AAri-values measured in e Tau's spectra observed 
before (filled tiangles) and after (open squares) T Tau. See 
text for details. 



B\it the a- value differs from zero by less than 2(T, which 
means that it is possible that the SIEs do not vary along 
the CCD. With this hypothesis we found the mean shifts 
are 11.8 ± 0.9, and 9.1 ± 0.7 mA for filled triangles and 
open squares respectively. The mean shift C derived for 
all points is 10.4 mA - dashed line in Fig.l - and standard 
error is ac =0-8 mA, or uc = 21 G, assuming .9 = 1 and 
Ao = 6400 A according to equation (1). The mean shifts 
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for filled triangles and open squares differ from C less than 
2(jc- So it appears resonable to use the mean C to remove 
SIEs, subtracting it from tlic; nK;asTu-(;d wavelength shifts 
between right and left circularly polarized spectra of T 
Tau. 

The corrected AAj.;-values in T Tau's spectra were 
used to calculate -By from equation (1). Resulting values 
of are presented in Table 2 for each of three observed 
spectra - the last three columns of the table. In some 
cases, a given line profile was distored due to cosmic ray 
events and thus was not used; these cases are marked with 
a dash. In fact, almost all absorption features in T Tau's 
spectra are blends of some number of lines, given in the 
first two columns. So we use an effective Lande factor, 

Qeff- 

derived from averaging over individual lines which form 
the blend. Here gi is the Lande factor of each line and di 
is its central depth. This approach appears resonable as 
long as the Zeeman splitting is less than the spectral res- 
olution of our instrument. We adopted T Tau's rf^- values 
from the VALD database specifing log g and Tg/ / of the 
star as input parameters. Averaging over all lines gives 
B\\ = +162 G with standard error crs = 30 G due to 
measurment uncertainty. Adding the error due to the un- 
certainty in C, the total uncertainty is: 

a ^ -^al + al ~ 37 G. 

Thus, we estimate B\\ = +160 ± 40 G from our 1996 ob- 
servations. 

3.2. 2002 observations 

The AAri (A)-valucs for all sets of e Tau's spectra are 
shown in Figure 2. It can be seen that there is no sta- 
tistically significant trend in these data around an aver- 
age value AXri = —0.02 ± 0.32 mA. Assuming geff = 1 
and Ao = 6700 A (the average wavelength of the spec- 
tral band), this value corresponds to B|| = —0.5 ± 8 G. 
Small values (AAw = -0.93 + 0.43 mA or = -20 + 10 
G) were also found for an additional giant star 75 Tau 
(K2III). Thus, we conclude that the adopted method of 
SIEs removal works well and confirmed a'posteriori our 
choice of e Tau as a zero-field standard during the 1996 
observations. 

To illustrate the reliability of our technique we in- 
cluded the magnetic star HD30466 in the target list of our 
2002 observations. Babcock (1958) found S|| = +2320 + 
340 G and B|| = +1890 + 130 G for HD30466 from two in- 
dependent observations. Profiles of the Cr i 6661.08 line in 
our HD30466's RCP and LCP spectra are shown in Fig.3: 
they are clearly shifted relative each other (AA,.( ~ 126 
mA), indicating ~ +2.0 kG, adopting g = 1.48 (Kupka 
et al. 1999) 

As can be seen from Table 3, relative shifts of lines 
between spectra with opposite polarizations in T Tau's 



Table 2. Results of 21.11.1996 observations of T Tau 
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spectrum are much less, and wc find — +140 ± 50 
G after averaging over all the data. The uncertainty is 
calculated by taking the standard error of the individual 
points. 

4. Discussion 

From observations in 1996 and 2002, using different meth- 
ods of SIE correction, we find similar (+160 ± 40 G and 
+140 + 50 G) values of the mean surface longitudinal mag- 
netic field in T Tau. As far as such small field values are 
close to the detection limit, it is very important to show 
that the measured wavelength shift is magnetic in origin. 
To test this we have measured relative shifts of the Sll 
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Table 3. Results of 22.02.2002 observations of T Tau 
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Fig. 2. AArj-values derived from three pairs of e Tau's 
spectra observed in 22.01.2002. See text for details. 
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Fig. 3. A part of RCP (thin line) and LCP (sohd hne) 
spectra of HD30466 
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6716.44 A and 6730.82 A forbidden emission lines in our 
2002 spectra, which prcsuinbly originate far from the star 
in a region with (nearly) zero magnetic field. We found 
AAr-i = +0.8 ± 1.4 mA for these lines which would give 
20 ± 35 G for .9e// = 1. This result can be considered as 
evidence that the measured wavelength shift of the stellar 
absorption lines is magnetic in origin. 

To further test if it is indeed possible to detect such 
low value from our spectra we performed the follow- 
ing test. The RCP and LCP spectra of ID stellar atmo- 
sphere were calculated with the SYNTHMAG program 
(Piskunov, 1998), assuming that \ogg and Te// are equal 
to that of T Tau and B = B\\ = 200 G, i.e. with magnetic 
field lines perpendecular to the stellar surface. Profiles of 
lines within 6300-6500 A spectral band were artifically 
broadened to take into account stellar rotation and spec- 
tral resolution and binned to the appropriate dispersion. 
Each of these artifical specra were used to produce three 
new ones by means of adding noise at the appropriate 
level. Finally we derived three pairs of RCP and LCP 
spectra with S/N-ratio near 100, as in the case of our 1996 
and/or 2002 observations. By means of cross-corellation of 
RCP and LCP profiles we then found B|| = 164 ± 32 G 
in agreement with the initial value of 200 G. Therefore, 
we belive that we indeed have recorded T Tau's surface 
magnetic field Sj| ^ 150 G, in spite of the relatively large 
scatter of the individual mesurements - see Table 2& 3. 



Regardless, our observations show T Tau's large scale B\\- 
value did not significantly exceed 200 G. 

This value is much less than the average value of T 
Tau's surface magnetic field (B> 2 kG) found by Guenther 
ct al. (1999) and Johns-KruU et al. (2000). This will re- 
sult if there are regions with opposite polarity of mag- 
netic field lines at the stellar surface. This could indicate 
the magnetic field of T Tau is non-dipolar. Note however, 
that even in the case of a magnetic dipole, the observed 
B[| = 0, if the dipole's axis is perpendecular to the line 
of sight. To demonstrate this more clearly we considered 
a magnetic dipole with equatorial field strength B= 1.6 
kG and calculated using the SYNTHMAG program the 
relative shifts of RCP and LCP profiles of the Fe i 6430.8 
line for different inclination angles j3 and then derived the 
"observed" longitudinal field. We found B|| ~ 900, 400 
and 200 G for /? 0°, 60° and 75° respectively. The mean 
unsigned surface field values for these inclinations are 2.4, 
2.2 and 2.1 kG respectively, i.e. close to values found by 
Guenther et al. (1999) and Johns-KruU et al. (2000). 

T Tau is believed to have a low inclination to its ro- 
tation axis based on ?;sin i (Hartmann & Stauffer, 1989) 
and rotation period measurements (Herbst et al., 1985). 
Therefore, one can expect that the magnetic axis is also 
close to the line of sight if it is aligned with or close to the 
rotation axis. Thus the low value of the average longitu- 
dinal magnetic field of T Tau we found can be interpreted 
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in two possible ways: 1) as an indication of the non-dipole 
character of the stellar magnetic field, or 2) as evidence of 
a large 30°) angle between the rotation and magnetic 
axes. 

5. Conclusion 

Weak circular polarization has been found in photospheric 
absorption lines of T Tau, indicating a moan longitudi- 
nal magnetic field ^ 150 G. While this value is near 
the detection limit of our apparatus, we belive that it is 
real. In any case one can conclude from our data that _B|| 
on T Tau does not significantly exceed 200 G, i.e. much 
less than surface magnetic field strength of the star (> 2.3 
kG) fomid by Guenther et al. (1999) and Johns-Krull et al. 
(2000). This could be the result of the non-dipole character 
of the magnetic field on T Tau or evidence of a large an- 
gle between the rotation and magnetic axes. More precise 
observations are necessary to distinguish between these 
possibilities. 
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